Stiffening Transition in Vicinal Surfaces with Adsorption 
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We study the vicinal surface with adsorption below the roughening temperature, using 
the restricted solid-on-solid model coupled with the Ising model. We calculate the step 
tension 7 and the step-interaction coefficient B by employing a variant of the density matrix 
algorithm. We find a "stiffening transition" at a temperature T s where B vanishes. At T s 
surface free energy has the form f(p) — /(0) = jp + (const. )p 5 + •• • (p: surface gradient), 
which differs from the well-known p-p 3 form. 

Below the roughening temperature Tr, the vicinal surface which is a slightly 
tilted surface relative to one of the facet plane of a crystal, is well described in 
terms of terraces, steps and kinks (TSK picture). Since systems whose "elemen- 
tary excitation" is an extended linear object (like the step), belong to the Gruber- 
Mullins-Pokrovsky-Talapov (GMPT) universality class^'i-i-i-i-i'S* we have the 
following form of the vicinal surface free energy (per projected area): 

/(p) = /(0)+7P + Bp 3 + 0(p 4 ), (1) 

where p (oc step density p) is the surface gradient. This p-p 3 form of expansion is 
characteristic of the GMPT universality class. Physically, 7 is the step tension, and 
B represents the effect of step-step interactions. To be precise, 7 and B depends 
on the mean running direction angle (relative to one of the crystal axes on the facet 
plane) of the steps, which we denote by 0; we should then write 7 = j(9, T) and 
B = B(9, T) (T: temperature). If we set up the xy coordinates on the facet plane so 
that the y-axis corresponds to 9 = (i.e., steps are running along the y-direction), 
the angle relates to the surface gradients p x (along the x-direction) and p y (along 
the y-direction) as 

Px = -ahP cos 9, p y = -a h p sin 9, (2) 

where ah is the height of a single step (for convenience, we adopt the unit where 
ah = 1 in the following). With (|2|), we can regard the free energy / as a function of 
the gradient vector p = (p x ,Py), allowing us to write / = f{p). For systems where 
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the step-step interaction is short-ranged, there exists the following universal relation 
between j(8, T) and BlB^Tl^^ (which leads to the universal Gaussian curvature 
jump at the facet edgeEP'cP't 



(3) 




n 2 (k B T) 2 
~ 6 7(0, T)' 

where k-g is the Boltzmann constant and j(9,T) defined by 

j(9,T)= 7 (0,T) + d 2 7 (6,T)/d9 2 



(4) 



is the step stiffness. How the long-ra nge (inverse-square) step-step interaction modi- 
fies the relation (^) is also known. Since the essential mechanism of the GMPT 
form ([]]) is the non-penetrability of steps, the universality of this form and also of 
the relation (|3|) are fairly "robust". However, for adsorbed vicinal surfaces, there may 
occur breakdown of the GMPT form at a temperature T s where B(9,T) vanishes, 
which we shall show in the present Letter. Since smallness of B is characteristic of 
low-temperature stiff surface, the anomalous behavior at T s may be called stiffening 
transition. 



face. 



is been 



own that adsorbed atoms often change properties of the sur- 
),§),@,0,[|),[§,|2§ 



Recently, the adsorption effect on the 
behavior of step has been s tud ied ex- 
perimentally.OS'B'BS'EI' In this 
Letter, to discuss the adsorption effect 
on the vicinal surface, we take the re- 
stricted solid-on-solid (RSOS) model© 
on the square lattice, coupled with the 
Ising spin system representing the ad- 
sorbed gas. In the RSOS model, we re- 
strict each nearest-neighbor (nn) height 
difference Ah to be Ah = 0, ±1, which is 
a reasonable simplification because con- 
figurations with large \Ah\ are energeti- 
cally unfavorable. 

We assume that the gas atoms are likely to adsorb at step edge positions^ 
and that the adsorbed atom modifies the ledge energy locally. The Ising spins are, 
then, located on the bonds of the square lattice where the RSOS model is denned; 
the Ising spins form a 45°-rotated square lattice. We further assume ferromagnetic 
interactions (attractive interaction, in the lattice-gas picture) with coupling constant 
J between the nearest-neighbor spins on the rotated square lattice. We assume 
simple linear modification of the ledge energy as e — > e(l — aa) (a: Ising spin). This 
modification leads to the interaction between the RSOS model and the Ising model. 
The Hamiltonian of the RSOS-Ising coupled model is, therefore, written as 



Fig. 1. (a) RSOS heights (h ly ...,h 4 ) and 
edge variables in the mapped vertex model, 
(b) Decorated- vertex-model representation 
of the RSOS-Ising coupled system. Ising 
spins are represented by circles. 



H= ^2 e ( 1 ~ aa b(i,j))\hi - hj\ - J a bO-b> 



(5) 
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where hi is the integer surface height at site i, e the "bare" ledge energy, <Jb{i,j) = ±1 
the Ising spin variable on the bond b(i,j) connecting the nn site pair < i,j >. We 
should note that the RSOS condition (\Ah\ < 1 for each nn site pair) is implicit in 

We analyze the model by the transfer-matrix method. For this purpose, we 
extend the well-known mapping between the RSOS model and the vertex model on 
the dual lattice,©'© to obtain a "decorated" vertex model (Fig.l). The decorated 
vertex model can again be regarded as a 6-state vertex model with 19 x 16 = 304 
non-zero vertex weights (304- vertex model). For approximate diagonalization of 
the transfer matrix._we employ the product- wavefunction renormalization group 
(PWFRG) method, EP'EZP which is a variant of the White's density matrix renor- 
malization group (DMRG) method ("infinite-system" algorithm, to be precise). 
The PWFRG is specially designed to obtain the fixed point (= thermodynamic limit 
of the systemEa- 1 ) of the DMRG efficiently. Fo r th e vicinal surface problem, we have 
verified the reliability of the PWFRG methodic ; even with a small number of "re- 
tained bases" (in the DMRG/PWFRG terminology), the method gives close-to-exact 
results. 

For calculation of j(9,T) and B(9,T) in ([!]), we take a similar approach as we 
took in Ref. 3C) where we obtained 7 and 7 for 2 x 1-reconstructed Si(00J_l . e!P> E^P 

We introduce Andreev field rjc3> along 
the x-direction to tilt the surface, 
by adding a term -r]Y,m,n( h (™>+ 1 ,n) ~ 
h( m , n )) m the Hamiltonian (||) ((m, n) is 
the position vector of the lattice site). 
In the vertex-model representation, the 
surface gradient p along the x-direction 
is just the thermal average of the verti- 
cal edge variable of the vertex model, 
which can easily be calculated from 
the fixed-point wavefunction obtained 
by the PWFRG. By sweeping the field 
77, we obtain a p — rj curve (actual calcu- 
lation is very similar to that of the mag- 
netization curve for spin chainsBB 
See also, Ref. 35)). From the standard argument using the Andreev's Legendre 
transformation p — > 77, f(p) — > f(r]) = f(p) — pq (77 = df(p)/dp), we have from (|l|) 
(with p = p) 

77 = 7 + 3Bp 2 + (higher order). (6) 

Hence, from the PWFRG calculation, we draw the 77 — p curve and perform the 
least-square fitting to obtain 7 and B. For actual fitting, we adopt the fitting form 



Fig. 2. PWFRG results for the p — (3rj curves 
(a = 0.5 and J = 0.15). Temperature of 
each curve is k B T/e = 0.3, 0.35, 0.4, 0.45, 
and 0.5, from right to left. 



77 = A + A 2 p 2 + A 3 p 3 + A 4P 4 



(7) 



since the coefficients ^3 and A4 may not be small. If the relation (|3|) holds, we obtain 
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the step stiffness 7(0) as 



7(0) 



7T 



2An 



(8) 



In Fig. 2, we show p-(3r] curves {(3 = 
1/ObT)) for a = 0.5 and J = 0.15 
calculated by the PWFRG, at some 
temperatures. In Fig. 3, we show tem- 
perature dependence of 7(0) obtained 
from (H). In contrast to the "normal" 
cases where 7 is a monotonically de- 
creasing function of temperature, our 7 
clearly exhibits diverging anomaly at T s 
(kBT s /e w 0.4). The step tension 7(0), 
on the other hand, behaves smoothly 
near T s , only with weak reentrance. At 
T s , the 7] — p curve is no longer fitted by 
the form ([?]); in stead, it is well fitted by 

77 = A) + A 4 p 4 + A 5 p 5 + ■■■, (9) 

to a non-GMPT 



which is equivalent 
form of expansion 



Fig. 3. Temperature dependence of step stiff- 
ness 7(0) (triangles) and step tension 7(0) 
(circles) (semi-log plot). Solid and broken 
lines are guide to the eyes. 



fip) = /(0) + IV + CsP 5 + (higher order). 



(10) 



We thus have found vanishing of B at a temperature T s in the RSOS-Ising coupled 
system, which we call stiffening transition. As for the p — rj curve, w^should have the 
change in the critical behavior, from the ordinary square-root typeK^ p ~ (i] — T] c ) 1 ^ 2 
(rj > T] c , 7] c = 7) to a new one p ~ (rj — rj c ) l / 4: . We should stress that the stiffening 
transition has nothing to do with the roughening transition, which takes place well 
above T s (we estimate k^T^/e ~ 1.35 in the present case). 

We propose two possible mechanisms of the stiffening transition: (1) "single- 
particle" mechanism, and (2) intrinsically many-body mechanism n The first one is 
based on the "universal free-fermion picture" of the vicinal surface,™ in a direct way. 
In the transfer-matrix treatment for the coarse-grained vicinal surface, the statistical- 
mechanical problem of many-step system is equivalent to finding the ground state 
of a free-fermion system in one dimension. In this view, the ordinary parabolic 
dispersion of the single-particle energy oj(k) ~ ak 2 (near k ~ 0) is the source of the 
/9 3 -term in ([l]) with B = o"/r 2 /3. Hence, the vanishing of B and appearance of the 
non-GMPT form ( |l0| ) can be simply interpreted as the change of u>(k) from fc 2 -type 
to /c 4 -type. Since 7 is inversely proportional to the (scaled) squared fluctuation width 
characterizing the step roughness,^ we should have step-smoothening behavior at 



The second mechanism is based on formation of "bound states" of steps. If we 
"integrate out" the Ising-degrees of freedom, we have an effective RSOS model with 
effective interactions. If this interaction amounts to attraction between the steps, 
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formation of "bound state" (or, "bound steps", "step bunch") may become possible. 
Suppose that steps are forming re-body bound states. We regard each bunch of steps 
as a "composite particle" and assume that the free-fermion picture itself is valid 
for the system of composite particles. By p n , *y n and 7„, we denote the density, 
formation free energy and the stiffness, of the n-body step bunch. Note that we 
have p = np n , and roughly 7„ ~ n.71 and j n ~ 7271. Ea' From a simple dimensional 
analysis, we see that in the expansion ([l]) of the vicinal-surface free energy f(p), the 
coefficient B scales as B ~ 1/n 4 . Hence, if n diverges then B vanishes; T s is the 
condensation temperature of the steps. 

To summarize, we have discussed the adsorption effect on the vicinal surface 
below the roughening temperature, in terms of the restricted solid-on-solid model 
coupled with the Ising model. By employing the pro duct- wavef unction renormaliza- 
tion group method, we have obtained temperature dependence of the step tension 
and step stiffness. We have found a stiffening transition where the step interaction 
coefficient vanishes, and at the same time, the vicinal surface free energy takes a 
different form from the Gruber-Mullins-Pokrovsky-Talapov-type form. 

We should note that anomalous increase of step stiffness has been observed ^ 
for boron (B)-doped Si(OOl). Similar phenomenon has also been^l^erwdEp for 
high temperature Si(lll) where Si adatom layer is known to exist. O'f^'ciP'EHP The 
stiffening transition we have found may have relevance to these anomalous behavior. 
In this Letter, we have restricted the analysis to special (but typical, in our view) set 
of parameter values of the model. To explore full parameter space is an important 
problem, which may lead to discovery of other interesting phenomena and may also 
lead to clarification of the mechanism of the stiffening transition. Models with other 
lattice structures, and/or with other types of couplings are also to be explored, for 
quantitative explanation of experiments. 
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